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Earth's analemma is the lopsided figure eight marked out over a year by the position of the Sun 
in the sky observed at the same clock time each day. It shows how solar time deviates from clock 
time. The analemma's shape results from the tilt of Earth's axis of rotation relative to the plane 
of its orbit around the Sun and from the elliptical shape of that orbit. This tutorial paper uses 
vector analysis of the Earth- Sun geometry and a numerically generated quantitative description of 
Earth's motion around the Sun to construct Earth's analemma. To visualize the geometry and the 
motion that give rise to the analemma is challenging, but this construction is a project within the 
capability of physics students who have had basic undergraduate mechanics. 



I. INTRODUCTION 




FIG. 1. Earth's analemma as a composite of pho- 
tographs of the Sun taken at the same location and 
time every day for a year. 

Gary Cooper's shootout in the movie "High 
Noon" almost surely was not at high noon. He 
was using clock time as established by the rail- 
roads. Even if he happened to be at the lon- 
gitude of his time zone corresponding to solar 
time, there are only four days in the year when 
the Sun crosses the meridian - its daily highest 
point in the sky - at 12:00:00 clock time. 

During a year solar noon varies from 16 min 
33 s earlier to 14 min 6 s later than clock time. 



This variation is apparent in Fig. [T] which is 
a composite of pictures taken in Cascade, Col- 
orado at the same clock time each day over the 
course of a year.^ The curve outlined by these 
365 images is Earth's analemma, and the pur- 
pose of this paper is to show how to construct it 
from basic principles of geometry and physics. 

The collection of images of the Sun in Fig. [l] 
is low in the western sky because the pictures 
were taken in the afternoon. The belt of the fig- 
ure eight, i.e., the points slightly below where 
the lines of the eight cross, moves along an arc 
that is the Earth's equator projected onto the 
sky, i.e., the celestial equator. You can imagine 
the figure eight of the analemma rising some- 
what tilted in the east - upper side first; then 
the figure eight, long axis perpendicular to the 
celestial equator, moves across the sky and sets 
in the west - upper side last. At noon the long 
axis of the figure eight lies along the meridian, 
half way between the eastern and western hori- 
zons. 

The analemma is a record of the Sun's posi- 
tion in the sky as observed from the same point 
on Earth at time intervals of exactly one day. 
The analemma displays the variation of the an- 
gle between the Sun's path in the sky and the 
celestial equator. Angle measured perpendicu- 
lar from the celestial equator is called declina- 
tion. The seasonal variation of the Sun's decli- 
nation from -23.44° to +23.44° and back again 
gives the analemma its long dimension. A small 
periodic variation in the time of day when the 
Sun crosses the meridian gives the analemma its 
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width of a few degrees P 

This paper shows how the relative motion of 
Sun and Earth results in the two angular varia- 
tions that produce Earth's analemma. The pa- 
per provides an analytic explanation in terms of 
two useful vectors, one that describes the mo- 
tion of the actual Sun, and the other that de- 
scribes a fictitious average or "mean Sun" that 
corresponds to clock time. These vectors are 
defined in section Ull 

Section [TTT| uses diagrams to show how geome- 
try results in the small periodic variation in the 
daily time at which the Sun crosses a merid- 
ian. These diagrams show how tilt produces an 
analemma even if the orbit is a perfect circle; 
they also show how to extract useful informa- 
tion from the two vectors. The vectors are then 
used to produce the analemma for the special 
case of a circular orbit. The result is a sym- 
metric figure-eight analemma with its mid-point 
exactly on the celestial equator. 



Section |IV| extends the vector approach to 
generate the analemma arising from an elliptical 
orbit. This section uses a numerical calculation 
of Earth's orbit and shows how its ellipticity 
gives the analemma the different sized loops and 
slightly lopsided shape you see in Fig[l] The ex- 
cellent agreement of the calculated analemma 
with the observed one is confirmed by laying 
one over the other. 

A higher level of vector analysis is exhibited 
in section [V| by using vector algebra to con- 
struct the analemma. The outcome is elegant 
and compact. 

In sections IVII and I VIII are a few fun facts 
for fostering facultative physics fora, some sug- 
gestions for student projects that could use ex- 
tensions of the ideas and techniques developed 
here, and some general conclusions. 



II. TWO VECTORS FOR 
CONSTRUCTING AN ANALEMMA 

It is the variation between clock time and Sun 
time that gives rise to the analemma. Conse- 
quently, to construct the analemma it is useful 
to define two vectors. One, rs(t), correspond- 



ing to the Sun's observed motion; the other, 
fm{t)^ corresponding to the motion of the fic- 
tional mean Sun and thus to the passage of clock 
time t. 

The vector rs(t) points at the real Sun as it 
makes its passage across the sky. The vector 
varies seasonally both in position and rate as 
it tracks the apparent geocentric motion of the 
Sun. 

The vector fm{t) points at the mean Sun. Be- 
cause the mean Sun defines the rate at which 
clocks run, rm(t) is in effect a clock vector. 
Thus, if frn{t) is set up on the equator at the lon- 
gitude used to define a time zone, it will point 
at the mean Sun's location directly overhead at 
exactly 12:00:00 noon in that time zone. For 
someone living in the U.S. Mountain Time Zone 
(UTC-07) that longitude is 105 °W. Set on the 
equator at that longitude, the vector fmit) will 
point directly overhead at the fictitious mean 
Sun as it crosses the meridian of people observ- 
ing from 105 °W longitude — for example, from 
downtown Denver, CO. 

As you will see, the analemma's shape is im- 
plicit in the difference between the two vectors: 



A{t) = rs{t)-rm{t). 



(1) 



Figure [2] shows the positions of the Sun and 
the mean Sun around a planet at successive 
times. (To make the tilt more evident, it has 
been drawn as 45°.) The dots are locations of 
the tips of the vectors fs{t) and fm{t) at 14-day 
intervals. Thus, the figure maps two spatial and 
one time dimensions into three spatial dimen- 
sions. The connected dots correspond to the 
angular position (two spatial dimensions) of the 
mean Sun exactly over a particular meridian — 
the meridian of the observer — day after day, 
and the arc length between dots of the mean 
Sun corresponds to an exactly integer number 
of days (the time dimension). The dots of the 
actual Sun are spaced by the same time interval 
as those of the mean Sun, but, as you will see, 
except in four instances the actual Sun can not 
cross the observer's meridian at the same time 
as the mean Sun. The dots of the real Sun move 
back and forth relative to the dots of the mean 
Sun. 
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FIG. 2. Intersecting circular paths traced out by 
the tip of the mean Sun (clock) vector — shown 
by connected dots — and the tip of the Sun vector 
for a circular orbit. The points of intersection of 
the two circles correspond to Earth-Sun equinoxes. 
The axes are oriented so the x-axis lies on the line 
of intersection (the line of nodes) of the two circles; 
the z-axis is perpendicular to the plane of the mean 
Sun and, therefore, parallel to the planet's axis of 
rotation. For visual impact the plane of rs is here 
tilted 45° to the equatorial plane swept out by rm; 
the dots are 14 days apart if the orbital period is a 
year. 



To extract the analemma from fs{t) and 
^m(0 you. need to represent them in a coordi- 
nate system. It is customary to use the pro- 
jection of Earth's system of spherical coordi- 
nates (latitude and longitude) onto the sky. For 
example, the circular arc formed by projecting 
Earth's equator onto the sky defines the celes- 
tial equator. 

Other projected lines of Earth's latitude de- 
fine corresponding lines of celestial latitude, i.e., 
declination. Lines of celestial longitude are de- 
fined similarly by projecting Earth's lines of lon- 
gitude onto the sky. 

Because of Earth's rotation, the lines of celes- 
tial longitude rise in the east, sweep overhead, 
and set in the west. By convention the line of 
zero celestial longitude lies in the plane defined 



by the Earth's axis and the line of nodes. 

To represent fs{t) and fm{t) it is convenient 
to use the celestial spherical polar coordinates 
in a Cartesian reference frame with its origin at 
Earth, its z-axis along Earth's axis of rotation, 
and its x- and y-diX.es in the equatorial plane. 
Referring to Fig. [2j you see that the x- and y- 
axes lie in the plane of the circle of connected 
dots traced out by rm(^)- 

By convention the x-axis is set to lie along 
the line of intersection (called the line of nodes) 
of the two circles and to point toward the point 
where the rising rs {t) of the tilted circle crosses 
the horizontal circle. In the geometry of Earth 
and Sun this crossing occurs around March 20 
— the spring equinox of the Northern Hemi- 
sphere. 

In this frame of reference fm{t) tracing out a 
circle in the x-y plane (Fig. [2| can be written 

{t) = COS ujtz -\- sin ujt j, (2) 

where z and j are unit vectors along the x- and 
y-axes respectively, t is time, and uj is the an- 
gular velocity of Earth on its axis. The time 
coordinate is measured from t = at the March 
equinox. 

The time dependence of fmit) is made to be 
such that after exactly one day, it again points 
directly overhead at the mean Sun as it crosses 
the observer's meridian. This time interval be- 
tween two successive occurrences of the mean 
Sun crossing any chosen meridian is defined to 
be 24 h, the internationally agreed upon clock 
length of an Earth day. Clocks are built to tick 
off exactly 24 h from one meridian crossing of 
the mean Sun to the next. 

It is important to understand that because 
Earth is traveling around the Sun, fmit) rotates 
through more than 27r radians in a day. This 
extra rotation occurs because even if Earth were 
not rotating on its axis, over the course of a year 
(365.25 d) the Sun and fm{t) would still make 
a rotation of 27r radians around Earth. 

To see the origin of this extra rotation, imag- 
ine Earth did not rotate on its axis at all. Then 
if on March 20 (for example) fmit) points at the 
Sun, six months later it will point directly away 
from the Sun, and after a year, it will once again 
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point at the Sun. Because of this extra rotation, 
in a year of 365.25 d an observer on Earth sees 
the Sun circle the Earth 366.25 times. Because 
of this extra rotation, the celestial longitude of 
the mean Sun observed from Earth advances by 
ujo = 27r/365.25 = 0.01720 rad each day. 

To take into account this extra rotation, 
fm{t) must rotate through an angle of 27r + 
27r/365.25 radian each day. Then after n days 
of rotation fm{t) will have rotated through an 
angle of n27r + ncjo- If the tip of frn{t) made a 
dot every 14 days, the dots would form a circle 
like the connected dots in Fig. [2] 

If Earth's orbit around the Sun were circular 
and if Earth's axis were not tilted relative to 
the plane of the orbit, the Sun would move in 
the sky exactly like the mean Sun; fs{t) would 
trace out a circle exactly like that traced out 
by r^{t). To see that this is so, relate fs{t) to 
fm{t) by expressing the time t as the number of 
days n. Then for no tilt and a circular orbit, 
the Sun vector is 

(^) = cos nujQ I + sin ncjo J- (3) 

The clock vector rotates faster, but, as discussed 
above, every 24 h it points at the mean Sun 
crossing the same meridian that it crossed a day 
earlier. It is the same meridian, but that merid- 
ian is aligned with an arc of celestial longitude 
cjo = .01720 rad ahead of the one it aligned with 
the day before. As a result, 

^m(^) = cos(n27r + ncjo) ^ 

+ sin(n27r + tiuoq) j. (4) 

Because cos(n27r+na) = cos(na) (and similarly 
for the sine), it follows that every 24 h rm(n) = 
fs(n), i.e., 

(n) = COS nujQ ^ + sin nuoQ j (5) 
rs(n). 

Thus, for an untilted circular orbit the 
analemma would be a single image of the Sun 
always at the same place in the sky at the same 
time each day. 



III. EFFECT OF AXIS TILT ON THE 
ANALEMMA 

Because Earth's axis of rotation is tilted 
23.44 ° relative to the plane of its orbit around 
the Sun, the equality of Eq. |5]does not hold for 
the real Sun. fs{t) traces out a path at an an- 
gle to the path traced out by fm{t). The tips of 
the two vectors trace out paths similar to those 
shown in Fig. |2] 

The analemma is shaped by the tilt of Earth's 
axis of rotation and by the ellipticity and ori- 
entation of its orbit. ^ The major effect comes 
from tilt. This is apparent when you construct 
an analemma for an orbit that is a circle. A 
tilted, circular orbit is also easier to deal with 
than an elliptical orbit, so it's a good procedure 
to first try out fs{t) and fm{t) in this simpler 
context and then modify them to take into ac- 
count the actual orientation and elliptical shape 
of Earth's orbit. 

Because of the tilt of Earth's axis, the Sun's 
angle relative to the celestial equator (its decli- 
nation) changes over the course of a year. Thus, 
for Earth, the points of fs{t) on the tilted cir- 
cle trace out the seasonal movement of the Sun 
as it goes from a declination of (5 = 0° at 
the March equinox to S = 23.44° (the Northern 
Hemisphere's summer solstice) to ° (autumnal 
equinox) to —23.44 ° (winter solstice) and so on. 
Figure [2] illustrates the variation in declination 
but with the tilt set to 45 °to make the effect 
more apparent. 

It is because its path is at an angle to the 
equatorial plane that the real Sun arrives over 
a given meridian at a different clock time each 
day. That is to say, when the mean Sun crosses 
a meridian, the real Sun will be over a slightly 
different meridian. Because Earth rotates at 
15° per hour, the angular difference between 
the longitudes of these two meridians corre- 
sponds to a time difference. For example, on a 
day when it is exactly clock noon at your merid- 
ian and the real Sun is overhead a meridian 1 ° 
to the east, it will cross your meridian 4 min 
later at 12:04 pm. A difference of 1 ° in longi- 
tude corresponds to 4 min of time, and 1 arc- 
min corresponds to 4 s of time. 
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FIG. 3. This is the appearance of the two intersect- 
ing circles shown in Fig. [2] as seen by a viewer look- 
ing directly down the z axis onto the plane of the 
mean Sun, the planet's equatorial plane — shown 
here by the connected dots. 
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FIG. 4. The mean Sun (o) moves at a steady rate 
around the celestial equator. The tilted orbit's pro- 
jection onto the plane of the mean Sun (Earth's 
equatorial plane) is shown by the +'s. The o's are 
the projection of the +'s onto the celestial equator, 
and they move back and forth relative to the mean 
Sun, i.e., relative to the regular time of a clock. 
(The tilt in this figure is 45 °.) 



Figures |3] and [4] show how the tilt causes the 
real Sun to vary in celestial longitude relative 
to the mean Sun. It is worth examining these 
figures closely both to understand the geometric 
origins of the effect and because they clarify how 
to extract a quantitative account of the effect 
from the two vectors fs{t) and fm{t)- 

This difference in the longitudes of the two 
meridians can be represented by a length of arc 
on the equatorial path of the mean Sun. To get 
this arc length, project fs{t) onto the equatorial 
plane, and then extend the projection to inter- 
sect the circle of the mean Sun. Figures [3] and 
|4] illustrate the construction of the arc length, 
and Fig. |4] shows how the movement of the Sun 
causes this arc length to vary as the real Sun 
moves back and forth relative to the meridian 
of the mean Sun. Combined with the seasonal 
change in declination this variation in arc length 
produces an analemma that is a figure eight. 

Why a figure eight? The geometry of Figs. [S] 
and |4] answers the question. These figures show 
the relation between the mean Sun and the ac- 
tual Sun over a succession of days. On the 
equinox an observer sees them both at the same 
meridian at the same time. But Fig. [3] shows 
that this alignment is soon lost. The figure is a 
view along the z axis perpendicular to the plane 
of the mean Sun, i.e., the equatorial plane. See 
how the tilt foreshortens the projection. This 
projection is shown by the +'s in Fig. |4j i.e., 
the +'s are the projection of the tip of fs onto 
the plane of fmit)- An observer will see an an- 
gular separation that corresponds to the arc be- 
tween the extension (the o's) of the projection 
of fs onto the arc along which the mean Sun is 
displayed (the o's). 

In Fig. [4] the time interval between successive 
o's, +'s, and o's is the same (~ 13 d), yet, as the 
figure shows, the o's move back and forth rel- 
ative to their corresponding o's.l^ To check this 
statement place a straight edge on the box (□) 
in Fig.[4]that marks the center of the graph, and 
notice that a line from the box to a + connects 
directly to a o. This o is the projection of the 
+ onto the equatorial circle. If you now follow 
the behavior point-by-point, you can see that at 
the equinox the o lines up with its correspond- 
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ing o — the mean Sun and the actual Sun cross 
the meridian at the same time. But then, at 
first, the Sun lags behind the mean Sun. These 
lags accumulate until the mean Sun is half way 
to the solstice (45 ° on the diagram); after that 
point the o moves faster than its correspond- 
ing o, and catches up to it so they both reach 
90 ° together. The o passes its o, but after 135 ° 
the o slows down and they both reach 180 ° to- 
gether. 

The cycle repeats over the second half of the 
year, so the back-and-forth motion occurs twice 
during the time that the up-down motion — the 
change in declination — goes through its cycle 
corresponding to one period of the orbit. Back- 
and-forth twice for up-and-down once makes a 
figure eight. 

The preceding argument and conclusions are 
implicit in Eq. [l] The declination angle S is 
the angle between A and its projection onto the 
equatorial plane. If /c is a unit vector in the 
direction of the Earth's axis of rotation, then 

cos (^^-5^ = k-fs 

sin (5 = /c • fs, (6) 

where k and fs are unit vectors. 

To extract a numerical result from Eq. |6j you 
need representations of fs and f^n in the same 
Cartesian coordinate system. For a frame with 
k parallel to Earth's rotation axis and z and j 
in the equatorial plane, frn is as given in Eq. |5j 
For rs, however, Eq. [3] is correct only for the 
frame designated by primes in Fig. [5j in which 
fs lies in the z-j plane which is perpendicular 
to the z' axis, and, therefore, perpendicular to 
k . 

The angle Q between k and k is the planet's 
tilt; Q = 23.44 °for Earth. For the case where 
one plane is tilted by rotating it an angle 
around the x-axis, the connection between the 
primed and unprimed unit vectors is, as you can 
see from Fig. [5) 
^/ 

z = z 

J = cos j -\- sin k 

k = — sinO j -\- cosO k. (7) 




FIG. 5. The connection between the unit vectors i, 
J, /c, and the primed unit vectors of a frame rotated 
an angle about their shared x-axis 

Using Eqs. [t] to express fs{t) = cos ooti -h 
sin ujtf in terms of the unprimed reference 
frame,!^ you will get 

fs{t) = cos UJtt 

+ sin ujt cosO j 

-\- sin ut sin 6 k . (8) 

Now you can evaluate Eq. [6] and find the dec- 
lination as a function of time: 

sin 6 = sin ut sin 
= sin nujQ sin 
S = sin~^ (sin nujo sin 0) . (9) 

Notice that, as expected, over the course of a 
year, i.e., from n = 1 to n = 365, 5 varies from 
+^ to —0 and back again. 

What about the variation of the difference 
along the equator? To find this follow the steps 
used to obtain Fig. |4j First, project fs(t) onto 
the plane of the mean Sun. The projection is 
the z and j components of Eq. [8) so you just 
omit the k component and get the projected 
vector fsxy where 

^Sxy(^) = COSUjtZ 

-\- sin ujt COS j (10) 
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This vector locates the +'s in Fig.[4j But you 
want to compare to the o's not the +'s but the 
o's, i.e., the points located by the extension of 
^Sxyi to the unit circle defined by f^. To make 
this extension, scale rsxy to make it be the same 
length as f^. Because f^n is of unit length, you 
scale fsxy by dividing by its magnitude |rsxy| 
\/ ^Sxy • ^Sxy , which from Eq. 10 



IS 



\^Sxy\ 



\/ cos^ ujt + sin^ ujt cos^ ^, (11) 



so that 



^Sx 



COS ujti-\- sin ut cos 6 j 
V cos^ ujt + sin^ ujt cos^ 



(12) 



Equation 12 describes a unit vector fsxy that 
is the same length as the unit vector f^. In 
Fig. [4] the tip of fsxy constructed in this way 
lies on the arc defined by fm{t)- 

The arc length on the unit circle equals (in 
radians) the angle, call it a, between and 
^Sxy You can find cos a from the scalar product 
of the two unit vectors 



ksxyl 

cos^ ujt + sin^ ujt cos 
\/ cos^ uot + sin^ uot cos^ 



(13) 



Alternatively, you can use the vector cross 
product to find sin a. 

Sina = k- {Vrn X ^Sxy) 

sincjtcosa;t(cos6> — 1) 



\/ cos^ ujt + sin^ ujt cos^ 9 



-sin(2cjt) sin^(-?^ 



\/ cos^ uot + sin^ uot cos^ 



(14) 



Using the cross product has some benefits. 
One is that when and fs are nearly equal, a 
is small enough so that sin a ^ a. Another is 
that application of the double-angle trig iden- 
tity makes it evident that the abscissa varies 
with twice the frequency of the ordinate. 

Equations |9] and [M] are parametric equa- 
tions of the analemma, and they produce the 
analemma curve shown in Fig. [6j This plot of 



Analemma from a tilted circular orbit 
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FIG. 6. Earth's analemma if Earth's orbit were a 
circle. Each dot is 1 d. 



declination 5 vs. a, the difference between the 
longitude of a Sun with a circular path and the 
longitude of the mean Sun, was made using Ex- 
cel. As foreseen, it is a symmetric figure eight. 
As is customary, the units of the abscissa are 
in minutes of time rather than in radians or 
degrees. 



IV. EFFECT OF ORBITAL 
ELLIPTICITY ON THE ANALEMMA 

Because Earth's orbit is not a circle but an 
ellipse with the Sun at one focus, the actual 
analemma (shown in Figs, [l] [IT] and 12) differs 
noticeably from the circular orbit's analemma 
shown in Fig. |6] The effects of the ellipticity 
can be included by replacing rs(t) for a tilted 
circular orbit (Eq. |8| with rs(t) for the actual 
elliptic orbit. 

A good way to get the fs that corresponds to 
the correct shape of Earth's orbit is to numeri- 
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cally integrate the equations of motion of Earth 
around the Sun to obtain rE(t), a vector that 
traces out Earth's orbit, and then use the fact 
that fs{t) = -Mt)' 




FIG. 7. Step 1: With the orbit (diamonds) in the 
plane of the mean Sun (connected dots) set the co- 
ordinate axes to have their origin at the focus of the 
ellipse and the y-axis along the ellipse's major axis 
2a. Perihelion is at y = — a(l — e). 

To do this integration you need the right ini- 
tial conditions. A good way to find them pro- 
ceeds in three steps. Step one, find the initial 
conditions in a coordinate frame aligned with 
the axes of the ellipse. This orientation of the 
axes makes determining the initial conditions 
particularly simple. 

Step two, transform the initial conditions de- 
termined in step one into the frame of reference 
that lines up its x-axis with the line of nodes 
— the line of intersection of the orbital plane 
with the plane of the circular orbit of the mean 
Sun (see Fig.[2|. For Earth, this corresponds to 
a rotation of the orbital plane 12.8° about the 
2;-axis. 

Step three, transform the initial conditions to 
correspond to the orbital plane tilted around its 
line of nodes. 

Numerical integration with the initial condi- 
tions from step three gives a vector te that ac- 
curately represents Earth's orbit. The resulting 




FIG. 8. Step 2: Rotate the axes to align the x-axis 
with the orbit's line of nodes. In this picture, the 
axes have been rotated around the z-axis by — 60 ° 
so the major axis of the ellipse is now at 60 ° to 
the ^-axis. The line of nodes is the segment of the 
X-axis crossing the ellipse. 




FIG. 9. Step 3: Rotate the axes around the x-axis 
so that the orbital plane is tilted 23.44 ° around the 
line of nodes (solid line). To make the tilt more 
apparent it has been set to 40 ° in this picture. 
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rs, which equals — te, is no longer the unit vec- 
tor of Eq. [3] sweeping out a circle. Now it sweeps 
out an ellipse with eccentricity e = 0.0167 tilted 
around the line of nodes at an angle of ^ = 
23.44° relative to the plane of the mean Sun, 
speeding up at perihelion and slowing down at 
aphelion. This vector fs, the mean Sun vec- 
tor Tm, and slightly modified versions of Eqs. |6] 
and [M] yield a quantitatively accurate represen- 
tation of Earth's analemma. 



-1.5 1 




FIG. 10. The sum of the distances from any point 
on the curve of an ellipse to the two points labeled 
"focus" is a constant equal to the length of the ma- 
jor axis 2a. The eccentricity e is the distance from 
the center of the ellipse to either focus, measured in 
units of the semi-major axis. 

Construction of this more elaborate fs vector 
requires some basic properties of ellipses sum- 
marized in Fig. 
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The two points labeled "fo- 
cus" are key. The curve is defined by the prop- 
erty that the distance from one focus to any 
point on the curve plus the distance from that 
point to the other focus is constant. (Use a ruler 
to check that this statement is true for Fig. 10 ) 
Another important feature is the eccentricity e. 
It measures the distance of the focus from the 
center of the ellipse in units of the semi-major 
axis a. (Use a ruler to confirm that the eccen- 
tricity of the ellipse in Fig. 10 is e = .75.) This 
value was chosen to be much larger than the 
0.0167 eccentricity of Earth's orbit to make the 
ellipticity more apparent in the figure. 

It follows from the general properties of an 



ellipse that 

semi-minor axis: b = a\/(^^-^^ (15) 
area of ellipse: A = irab. (16) 

Newton's law of universal gravitation states 
that the force F21 exerted on a spherical mass 
m2 by a spherical mass M is 



F21 



GMm2 (^2 — fi) 
1^2 - rip |r2 - ri| 



(17) 



where fi and r2 are the position vectors of the 
two masses, and G is the constant of universal 
gravitation. This equation is the inverse square 
law of gravitational attraction multiplied by the 
unit vector (r2 — ri)/(\r2 — ri|), which points 
along the line between the centers of the two 
masses. The unit vector asserts that this is a 
central force, and the minus sign at the begin- 
ning of the equation points the unit vector from 
m2 to M to refiect the fact that M is attracting 
m2. 

For M » 7X12 and with r = r2 — fi, Newton's 



laws of motion and Eq. [T7| lead to a differential 
equation for the orbit: 



-CM 



\ff \r\' 

The orbiting mass m2 has energy E where 

E 1^ ^ GM 
1712 2 \r\ 



(18) 



(19) 



and angular momentum L such that 



L 

1712 



r X V = r uj, 



(20) 



where Co is the angular velocity of the planet 
around the Sun. E and L are conserved (con- 
stant) quantities. 

For initial values of position and velocity such 
that the energy E < 0, the solutions to the 
differential equation are orbits that are ellipses 
with semi-major axis a and the Sun at one fo- 
cus. The following equations connect the ge- 
ometry of the elliptical orbit to the underlying 
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TABLE I. Some Parameters of Earth's Orbit 



e 


Eccentricity 


0.0167 


a 


Semi-major axis 


1 AU 


T 


Period 


1 y 


E 


Energy per solar mass 


-19.739 AUV 


L 


Angular momentum 
per solar mass 


6.2823 AUVy 





Angle from line of 
nodes to perihelion 


-78.2 deg 



physics: 



GM 



2E 
GM 



a 



(21) 
(22) 
(23) 



where T is the period of the planet's orbit. ^ No- 
tice that with a measured in AU and T in years, 
Eq. [23] becomes = 1 = 47r VGM for Earth, 
so that in these units GM = 47r^ when M is the 
mass of the Sun. 

These equations are interesting for what they 
say about planetary motion. They are essential 
for checking the validity of orbits produced by 
numerical integrations. They are also useful for 
obtaining initial conditions for orbits with pa- 
rameters different from those of Earth's orbit. 

To find initial conditions that correspond to 
Earth's orbit, that is, to an orbit that has e = 
.0167 and the other parameters given in Table [T| 
start with step one of the procedure outlined 
above. 

In a Cartesian coordinate frame on the ellipse 
with the origin at one focus, the x-axis coinci- 
dent with the ellipse's minor axis, and the 7/-axis 
with its major axis. Earth's position vector at 
perihelion is te = £ — a(l — e) j and its velocity 
is = Vpi^Oj. 

The quantity Vp is the magnitude of Earth's 
velocity at perihelion. You can determine 
using Kepler's law of equal areas in equal times 
(conservation of angular momentum L). 



A vector f from the Sun to a point moving 
with velocity v sweeps out area at a rate 



dA U ^ 

— — = -r X V. 
dt 2 



(24) 



This is essentially the formula for the area of 
a triangle — 1/2 the base times the altitude 
where the base is the arc length swept out in a 
time dt. Kepler and Newton tell us that for a 
central force, the quantity dA/dt is constant in 
both magnitude and direction. It follows that 
for Earth the value for that constant must be 
the same as the total area swept out in a year. 



i.e., Eq. 16, divided by the time of one year: 



dA 



dt 



nab 

~T~ T 
3.141 AUVy. 



(25) 



But, as Eq. 20 states, f x v \s also the angular 
momentum L of a unit mass. (For any given 
initial conditions, the orbit is independent of 
the mass of the orbiting object, so set m2 = 1 
to simplify the discussion.) 



L 



L 

^Earth 



2— 

dt 

6.282 AUVy 



(26) 



And, because at perihelion te and ve are ex- 
actly perpendicular. 



^Earth Ci{l - e)\vE\ 

= (1 - .0167) 

= .9833 Vp = 6.282, 



(27) 



and the magnitude of the velocity at perihelion 
is 



Vp = 6.389 AU/y. 



(28) 



In the chosen reference frame the position of 
Earth at perihelion is te = — .9833 j AU, and 
its velocity is = 6.389 z AU/y. It would be 
equally convenient to determine the initial con- 
ditions at aphelion. 

With these initial conditions an Euler integra- 
tion of the differential equation with an initial 
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half-step gives a good result. This is the method 
that Feynman used in his lectures."^ In the units 
used here a convenient time step for the inte- 
gration is At = 1/365.25 = .002738 y. This is 
the duration of 1 d (one day) so 366 steps of 
integration produce a complete orbit, and with 
this time step you can use the representation 
of Tm given in Eq. [5] The integration can be 
done with any mathematical software package. 
A spreadsheet works well, and the calculation 
used here is available with the online version of 
this paper. 

TABLE II. Initial values of te in AU and Vp in 
AU/y for calculating Earth's orbit in three different 
orientations 



Orientation of el- 
lipse relative to 
coordinate axes 


Perihelion on — ^ axis 


Orbit plane rotated -12.8° 


Orbit plane tilted 23.44° 







0.2184 


0.2184 


yo 


-0.9833 


-0.9588 


-0.8796 










0.3814 


VxOl/2 


6.3890 


6.2171 


6.2171 


VyOl/2 


0.05589 


1.4733 


1.3517 


VzOl/2 








-0.5643 



Integration with these initial conditions starts 
from perihelion (January 3) and gives the day- 
by-day position of Earth on its orbit with the 
specified eccentricity. However, this orbit needs 
to be correctly oriented relative to the idealized 
orbit of the mean Sun. 

Do this with steps two and three. Step two: 
make the line of nodes coincide with the x-axis 
by rotating the coordinate frame about its z- 
axis until Earth crosses the x-axis on March 20, 
the vernal equinox of the Northern Hemisphere. 

To find the angle of rotation use orbit data 
produced by numerical integration with a tilt 



TABLE III. Selected values of te from the numer- 
ical integration of Earth's orbit. 



i^cLie 




Vn 


Z"n 




3- Jan 


0.2184 


-0.8796 


0.3812 


0.9833 


1-Feb 


0.6638 


-0.6682 


0.2897 


0.9854 


20-Mar 


0.9959 


-0.00647 


0.00280 


0.9960 


21-Mar 


0.9962 


0.00940 


-0.00407 


0.9963 


20-Jun 


0.02009 


0.9323 


-0.4042 


1.0163 


21-Jun 


0.00317 


0.9325 


-0.4043 


1.0164 


22-Jun 


-0.0138 


0.9325 


-0.4043 


1.0165 


5-Jul 


-0.2317 


0.9084 


-0.3953 


1.01685 


22-Sep 


-1.00366 


0.01344 


-0.00583 


1.00377 


23-Sep 


-1.00349 


-0.00229 


0.00099 


1.00349 



angle but no rotation of the orbit plane. These 
data show Earth crossing the x-axis 13 d later 
than the actual March 20 equinox. This re- 
sult means you must rotate the reference frame 
clockwise by -13/365.25 x 360 = -12.8° rel- 
ative to the orbit to place the x-axis so that 
Earth crosses it on March 20.^ 

To rotate the coordinate frame, use the fol- 
lowing version of Eq. [7] modified to take into 
account that this rotation is about the z-axis!^ 

Xo = Xq cos a + ^0 ^ 

^0 = + (-0.9833 X -0.2221) = 0.2184 

7/0 = —^0 sin a + ^0 ^ 

^0 = + (-0.9833 X 0.9750) = -0.9588. 

These equations produce the second set of ini- 
tial conditions in Table [Hi 

Step three: tilt the orbital plane. Do this by 
applying Eq. [7|to the set of orbital coordinates 
labeled "rotated" in Table [n| you will get the 
initial conditions labeled "tilted" in Table |TT1 
Integration with these values gives Earth's orbit 
tilted relative to the equatorial plane. 

The result of the numerical integration is a 
table of values of x^, Vn^ and z^, the Cartesian 
components of rE(ncJo), the vector that locates 
Earth relative to the Sun on day n since peri- 
helion: 

rE(ncJo) = Xni^ynJ^Vnk (29) 
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Some results from the numerical integration 



are given in Table III These should make you 
confident that the calculation is accurate. Be- 
tween March 20 and 21 yn and Zn pass through 
zero as they should at the vernal equinox of the 
Northern Hemisphere. And between June 21 
and June 22, when the summer solstice occurs, 
Xn changes sign as it should. The entry dated 
July 5 was selected because on that date the 
distance has its largest value; in other words 
the calculation gives the correct date for Earth's 
aphelion. 

But notice that the aphelion distance is 
slightly larger than a(l + e); this discrepancy 
reflects the limitations of the method of numer- 
ical integration and the step size used here. It's 
a warning that integrations over more than a 
year may need a smaller step size or a more ac- 
curate method of integration. 

To construct the analemma, you need fs the 
vector that points from Earth to Sun. It is 

fs{nuJo) = -fE{nuJo) 

= -Xn l-ynj-ynk (30) 



To extract a and (5, modify Eqs. [TT] [T2j and 
[Mjto be Eqs.|3lJ|32l and[34]and use them to get 
a and 5 as functions of time, i.e., as functions 
of n the day number. 



\rsxy{nuJo) \ = ^xl^yl, 



so that 



Then 



so that 



fsxy(na;o) 



Xji z -\- y^i J 



sin a = k-frn X ^Sxy 



. -1 / —yn COS ncjo + sm nuo 
an = sm — 



(31) 
(32) 

(33) 

(34) 




FIG. 11. Here the calculated analemma is displayed 
overlaid on the observed analemma of Fig. [l] For 
the best overlap the calculated analemma has been 
rotated through 35.6°. This corresponds to about 
2:15 pm MST in reasonable agreement with the 2:28 
pm MST when the photographs of the observed 
analemma were taken. 



and 



Equation 34 is not entirely satisfactory for ex- 
tracting a. It is set up so that when n = 0, 
fs points 12.8° counterclockwise from the nega- 
tive y-axis while is pointing along the x-axis. 
This makes the angle between them 90 — 12.8 = 
77.2°. In principle this large angle is not impor- 
tant; the variation in a is still the correct effect; 
it's as though you were measuring the time at 
which the Sun is near its zenith with a clock 
in a time zone some seven or eight hours away. 
However, it is tidier to set f^n so that it corre- 
sponds to the direction of fs. To do this, add a 
phase constant (j) = —77.2° so frn becomes 

fm{n) = cos(nco'o + (/>) ^ + sin(nco'o + </>) J (36) 

and 

Xn sin(na;o + 0) - 2/n cos(ncjo + 0) 
^ . [o7) 



Sn = sin ^ 



(^35 j The analemmas shown in Figs. pT] an d are 

plots of S vs. a obtained from Eqs. |35| and 
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Earth's analemma 



-20 -15 -10/ -5 



deviation from mean 
solar time (mins of 



5 ••, 10 15 20 25 




FIG. 12. This is the same curve as in Fig. |11| but 
oriented to be at noon and with the scale of the 
abscissa expanded by x6. Each dot is a different 
day. 



37 for 



1 to n = 365 with 



-77.2° 



In Fig. 11 the calculated analemma has been 
laid over the observed analemma; the agree- 
ment between the two is quite satisfactory. The 



analemma in Fig. [12] is the same as in Fig. [TT] 
but with the horizontal dimension enlarged by 
6x. It shows the day- by-day variation of the 
arrival of the noontime Sun. 



V. VECTORIA'S SECRET 

The construction of the analemma in the pre- 
ceding sections emphasizes visualization of how 
the analemma results from the relative motion 
of Earth and Sun in three dimensions. The goal 
was to show how the geometry connects with 
the physics. 

But more formal, less visual approaches are 
possible, and they have their merits. Done 



properly they can carry you past various diffi- 
culties of establishing relative phases, of getting 
the algebraic signs correct, and they can pro- 
vide an aid or an alternative for those who find 
it difficult to visualize three dimensions. The 
following paragraphs make a fuller use of vec- 
tor algebra to reprise the construction of the 
analemma. 

The key question in constructing the 
analemma is "What is the angle a between fm 
and the projection of fs onto the plane of rm?" 
One way to find the answer is to use the fact 
that a is the dihedral angle between two planes, 
one defined by k and and the other by k and 
fs. Figure 13 shows the configuration. 




FIG. 13. The vectors k, rm, and fs define two 
planes. The variation between the positions of the 
Sun and the idealized mean Sun causes a to vary 
from day to day. 



The angle a between the two planes is the 
same as the angle between unit vectors drawn 
perpendicular to the two planes. This fact sug- 
gests the following procedure to find sin a: Con- 
struct unit vectors hs and hm normal to the two 
planes; take their cross product, which will give 
a vector in the k direction with magnitude of 
sin a; take the dot product of k with the cross 
product to find sin a (see Eq. 40). 



14 



Use cross products to construct hs 
k X fs 



ns 



and hn 



I A; X rsl 



k X Tyh 

\kxrm\ 
/c X r m 



(38) 



(39) 



where Eq. 39 follows from the fact that | k x 
Tml = 1 because k and fm are both unit mag- 
nitudes and perpendicular to each other. 

This is enough information to determine sin a 
by direct substitution into 

sin a = k ■ {hm x ns), 



(40) 



but using the vector identity A-{B xC) = {Ax 
B) • C and the fact that rm is 



yields the more compact expression: 

sin a = k • {n^ x hs) 
= {kx hm) • ns 



(41) 



(42) 



Eq. 37 



Substituting Eqs. 36 and 38 into Eq. [42] gives 



The vector expression for the declination was 
already used to obtain the result given in Eq. 35 



VI. FACTOIDS AND PROJECTS 

The analemma occurs because solar days vary 
in length. A curious consequence is that the 
longest solar day, about 24 h 30 s, occurs just 
after the winter solstice, the day on which in 
the Northern Hemisphere there are the fewest 
hours of daylight. Thus, "the shortest day of 
the year is also the longest day of the year" — 
or nearly. It is also interesting that for about 
three weeks after the solstice sunrise continues 



to get later. Nevertheless, the amount of day- 
light increases in the Northern Hemisphere be- 
cause sunset is getting later by an amount that 
more than offsets the later sunrise. The New 
York Times finds these factoids fit to print. 

This tutorial can be the basis for further 
projects. A project might be as modest as 
an exercise in vector manipulation. For exam- 
ple, you might replace hm and fs in Eq. 40 
with their cross product definitions from Eqs. 38 
and 39 and use the identity A x (Sx C) = 
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B{A-C) -C {A- B) to derive Eq. 

Or you could measure the declination of the 
Sun and observe the time at which it reaches its 
zenith and see if these agree with the predictions 
of your analemma. You will need to understand 
how to make f^n correspond to your clock time 

- a worthwhile task. 

You can modify the input parameters for the 
numerical calculation of Earth's orbit and ex- 
amine how Earth's analemma will change as the 
equinoxes precess. You can also replace Earth's 
orbital and rotational parameters with those of 
Mars and find the analemma that an observer 
on Mars will see. By properly modifying the 
input parameters of the numerical integration, 
you can find the analemma of any planet. 

The analemma is the result of the superposi- 
tion of two periodic motions at right angles to 
each other. In effect it is an interference pat- 
tern, a Lissajou figure. The rather small eccen- 
tricity of Earth has a surprisingly large effect on 
the analemma. Can you invert the argument, 
and find the eccentricity of Earth's orbit from 
its observed analemma? 

A nice project would be to calculate the tran- 
sits of Venus using this paper's approach. Nu- 
merically calculate the position vector fy of 
Venus relative to the Sun. Then construct a 
pointer from Earth to Venus tev = Ty — ve- 
Look at the angular difference between fs (= 
— te) and fEv as a function of time. Find the 
times at which the angular difference is between 

— 1/4° and +1/4°, i.e., approximately within 
the angular diameter of the Sun. You will need 
a method of numerical integration that is sta- 
ble over the long interval of time between the 
occurrence of one pair of transits and the next 
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pair. The Euler method wih not work, but a 
fourth order Runge Kutta method available on 
most packages of mathematics software may do 
the job. 



dents to go beyond praise for the achievements 
of Newtonian physics. This tutorial enables 
them to recreate and experience some of that 
achievement. 



VII. CONCLUSION 
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^ Wojtek Rychlik , (2005) , |< http : //epod . usra . ' 
edu/blog/2005/Ol/colorado- analemma . html> 

shows a composite of over 300 images taken at 
14:28:00 local time in Cascade, CO almost every 
day from Dec. 4, 2003 to Dec. 4, 2004. 

^ Astronomers call this small variation the "equa- 
tion of time," using the word "equation" in 
an obsolete sense meaning " correction." The 
analemma provides the correction to convert 
sundial time to clock time and vice versa. 

^ Astronomers call the tilt angle the "obliquity" of 
Earth's rotation axis. 

^ The number of o's is the same as the number of 
+'s is the same as the number of o's. 

^ You can write down Eq. [s] directly without go- 
ing through the transformation, but it is good 
to use this transformation here because it will 
be used two more times in constructing Earth's 
analemma. 

^ The conversion of a from radians to minutes of 

, . . 1 180 deg 60 min/h ^ 

time IS 1 rad — tt — 229.2 mm. 

TT rad 15 deg/h 

^ Keith R. Symon, Mechanics (Addison Wesley, 

Reading, MA, 1953) pp. xiii, 358 (see pp. 111- 



117). 

Richard P. Feynman, Robert B. Leighton, and 
Matthew Sands, Feynman Lectures on Physics^ 
Vol. 1 (Addison Wesley, Reading, MA, 1963) pp. 
9-6 - 9-9 

You can also use the fact that perihelion for 
Earth does not occur on the day of the sol- 
stice, December 21, as would be the case if 
you did not rotate the axes. Perihelion oc- 
curs 13 days later on January 3, so the ro- 
tation of the axes that puts perihelion at the 
proper date is (as it should be) the same as 
the rotation that places the March equinox 
on the X-axis. Another alternative is to look 
up the celestial coordinates of perihelion; see, 
for example, <http://nssdc.gsf c.nasa.govTl 
planetary/f act sheet /earthf act . html> 
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from Eq.[7|by permut- 



You can generate Eqs. 
ing X, 2/, and z (and j, and k. 
Joseph S. Chalmers, "Ground trajectories of 
geosynchronous satellites and the analemma," 
American Journal of Physics 55, 548-552 (1987). 



